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ABSTRACT 



o r We have studied the location of narrow-line Seyfert 1 (NLS1) galaxies and 

broad-line Seyfert 1 (BLS1) galaxies on the Mbh — cr relation of non-active galax- 
+3 \ ies. We find that NLS1 galaxies as a class - as well as the BLS1 galaxies of our 

comparison sample - do follow the Mbh — cr relation of non-active galaxies if we 
use the width of the [SII]6716,6731 emission lines as surrogate for stellar velocity 
dispersion, cr*. We also find that the width of [OIII]5007 is a good surrogate for 
a*, but only after (a) removal of asymmetric blue wings, and, more important, 
after (b) excluding core [OIII] lines with strong blueshifts (i.e., excluding galaxies 
which have their [OIII] velocity fields dominated by radial motions, presumably 
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outflows). The same galaxies which are extreme outliers in [OIII] still follow the 
Mbh — o~ relation in [SII]. We confirm previous findings that NLS1 galaxies are 
systematically off-set from the Mbh — cr relation if the full [OIII] profile is used to 
measure a. We systematically investigate the influence of several parameters on 
the NSL1 galaxies' location on the Mbh — cr plane: [OIII] cor e blueshift, L/LEdd, 
intensity ratio FeII/H/3, NLR density, and absolute magnitude. Implications for 
NLS1 models and for their evolution along the Mbh — cr relation are discussed. 

Subject headings: galaxies: active - galaxies: bulges - galaxies: evolution - 
quasars: emission lines 



1. Introduction 

Correlations between the properties of galaxies and the masses of their central black 
holes are of fundamental importance for our understanding of galaxy formation and evolu- 
tion. The relations between black hole mass and bulge stellar velocity dispersion, Mbh — 
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(Ferrarese & Merritt 2000, Gebhardt et al. 2000, Merritt & Ferrarese 2001) and between 
black hole mass and bulge mass (Marconi & Hunt 2003, Haring & Rix 2004) imply a close 
link between black hole and bulge formation and growth (e.g., Haehnelt & Kauffmann 2000). 

Mathur et al. (2001) emphasized the importance of studying the locus of NLS1 galaxies 
on the M B h— cr plane. NLSls are a subclass of Active Galactic Nuclei (AGN) characterized by 
small widths of their broad Balmer lines (e.g., Veron-Cetty et al. 2001). They are believed 
to have low-mass black holes (BHs) accreting close to the Eddington limit (e.g., Boroson 
2002). As such, they are an 'extreme' class of AGN in the local universe, potentially holding 
important clues as to the origin of the Mbh — cr relation (see Sect. 1 of Li et al. 2006 for 
a review on models), and whether it holds at all times and for all types of galaxies, and at 
low BH masses. 

However, previous studies of NLS1 galaxies, employing different samples and methods, 
led to partially conflicting results regarding the galaxies' location on the Mbh — cr relation. 
When using the width of the narrow [OIII]5007 emission line, FWHM[ ni], as replacement for 
stellar velocity dispersion a*, most authors found that, on average, NLS1 galaxies are off the 
M BU -o- [om] relation (Mathur et al. 2001, Grupe & Mathur 2004, Mathur & Grupe 2005a,b, 
Bian et al. 2006, Watson et al. 2007), while in the study of Wang & Lu (2001) most NLSls 
are, on average, on the relation (with a significant scatter in all cases). Furthermore, NLS1 
galaxies were found to be off (Wandel 2002, Ryan et al. 2007) or on (Botte et al. 2004) the 
black-hole-mass - bulge-luminosity, Mbh - -^buige, relation; and on (Botte et al. 2005, Barth et 
al. 2005) or off (Zhou et al. 2006) the Mbh — °~* relation. Most larger samples are actually 
based on FWHM[om]/2.35 as surrogate for a* (Nelson 2000), because the bulge velocity 
dispersion cr* is very difficult to measure in NLS1 galaxies due to superposed emission-line 
complexes. 

Given the potentially far-reaching consequences of a deviation of low BH mass, highly 
accreting AGN from the Mbh — cr relation, independent studies of large (r) samples of NLS1 
galaxies are important. No matter whether they turn out to be on or off the Mbh — cr relation 
of non-active galaxies and BLS1 galaxies, any result has important consequences: if NLSls 
lie off the relation, for the formation of the M B h — cr relation (how do objects move onto 
the relation ?) and/or the gaseous kinematics of their NLRs. If they lie on the relation the 
question is raised how to preserve Mbh — cr along the relation in objects which are rapidly 
growing their BHs but are typically in non-merging galaxies (Krongold et al. 2001). 

We have used an independent, large, homogeneously analyzed sample of NLS1 galaxies 
(55 objects) plus a comparison sample of BLS1 galaxies (39 objects), in order to study 
their locus on the Mbh ~ cr plane. We conclude that the width of the [SII] emission line, 
and of the narrow core of the [OIII] line (but only after removal of [OIII] lines with high 
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blueshifts), are the best available surrogates for stellar velocity dispersion cr*. Using these 
lines puts NLS1 galaxies (and BLS1 galaxies), as a class, on the M B h — cr relation of non- 
active galaxies. Throughout this paper, we use a cosmology with H =70 kms" 1 Mpc -1 , 
f2 M =0.3 and fl\=0.7. 

2. Sample selection and data analysis 

We use the sample of NLS1 and BLS1 galaxies of Xu et al. (2007). Sample selection, 
data preparation, and data analysis methods are described in detail in that work. In brief: 
the sample consists of NLS1 galaxies from the catalogue of Veron-Cetty & Veron (2003), 
plus a comparison sample of BLS1 galaxies from Boroson (2003) at z < 0.3, which have 
detectable (S/N>5) [SII] emission in their Sloan Digital Sky Survey (SDSS) DR3 (Abazajian 
et al. 2005) spectra. The galaxies of the two samples have similar redshift and absolute 
magnitude distributions. Xu et al. (2007) corrected the SDSS spectra for Galactic extinction, 
decomposed the continuum into stellar and AGN components, and then subtracted the stellar 
continuum contribution and the Fell complexes from the spectra. Emission lines were fit 
with Gaussian profiles. 

Special attention was paid to [OIII] profile shapes. The deviation of [OIII] profiles from 
single Gaussians is well known, and NLS1 galaxies show strong blue wings particularly often 
(Boroson 2005 and ref. therein). Authors who previously used qoni] from NLS1 galaxies were 
aware of this complication, and most of them did correct for it. In most cases they did so by 
using only the 'red' half of the [OIII] line profile to measure the FWHM. Wang & Lu (2001) 
decomposed the [OIII] profiles into core and wing. We followed two approaches in measuring 
FWHM[oni]: (a) determination of the total FWHM of the summed two-component Gaussian 
from a careful decomposition into a core and a blue wing component whenever present, 
FWHM [0 iii] totl ; and (b) the FWHM of the core component of [OIII], FWHM [0 ni] core . Emission 
lines other than [OIII] are well represented by a single Gaussian. The lines [SII]6716,6731 
(hereafter referred to as [SB]) were both fit with the same width. We follow the standard 
procedure of using the radius (i?BLR,)Tuminosity(L) relation of BLS1 galaxies (which appears 
to hold for NLS1 galaxies; Peterson et al. 2000) to compute BH masses, following Kaspi 
et al. (2005). We then compare our results with the Mbh — o"* relation of Tremaine et al. 
(2002; their eqn. 1 and 19) and of Ferrarese & Ford (2005, FF05 hereafter; their eqn. 20). 
An error in black hole mass of 0.5 dex (as in Grupe & Mathur 2004) was conservatively 
assumed. Formal measurement errors in FWHM are typically 5-20%. 

We note that all SDSS observations of our sample include emission from the entire NLR 
(so we are not systematically excluding the outer parts of the NLR in the nearest objects). 
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At the lowest redshift in our sample, z=0.034, an SDSS fiber diameter of 3" corresponds to 
2 kpc, larger than typical NLRs. 

3. Results 

We investigated how the following factors affect the galaxies' location on the M BH — cr 
plane: (a) different methods of fitting the [OIII] profile; (b) the use of other emission lines as 
surrogates for cr*; and (c) the influence of other parameters, like [OIII] blueshift and L/L E dd- 

3.1. NLS1 galaxies on the M BH — o"[om] plane 

First, we confirm that NLS1 galaxies as a class deviate from the M BH — cr relation of 
non-active galaxies if we use the full profile of [OIII] to measure cr (Fig. 1). We computed the 
deviation Acr := \oga obs — \oga pred of the NLS1 and BLS1 galaxies from the FF05 relation. 
Applying a Kolmogorov-Smirnov (K-S) test shows that the two samples are drawn from 
different parent distributions (K-S probability of being drawn from the same population: 
0.002). NLS1 galaxies follow the M BB — cr relation much better if only the narrow core of 
[OIII] is used to measure cr (Fig. 1). However, there are still a number of NLS1 galaxies 
with very broad [OIII] profiles and these lie systematically below the M BB — cr relation. In 
Sect. 3.3 we identify the cause of this. We show that the majority of these '[OIII] outliers' 
are characterized by high NLR outflow velocities, and therefore, their velocity fields are not 
dominated by the bulge potential and their line widths cannot be used as surrogates for cr*. 

3.2. Widths of other NLR lines as surrogate for cr* 

The strongest (high-ionization) emission line in AGN spectra is [OIII] 5007, which is 
commonly used for cr estimates. Because of its strength, we are able to decompose the line 
into a core and blue wing. Low-ionization emission lines generally do not show this spectral 
complexity, and we therefore also consider using their widths as surrogates for cr*. The 
strongest lines in the commonly available wavelength range between ~ 4000 — 7000A are 
[NII]6584,6548, [SII] and [OI]6300. [Nil] is, however, blended with Ha except in high- 
resolution spectra (see Zhou et al. 2006 for the use of (7[nii] as a surrogate for cr*). We focus 
on [SII] here, which is detected in all galaxies of our sample. As shown by Greene & Ho 
(2005), <7[sn] is indeed a good substitute for cr* in Seyfert 2 galaxies. 

Using the width of [SII] as a surrogate for cr*, NLS1 galaxies as a class are now on the 
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.Mbh — 0" relation (Fig. 1). Repeating the K-S test on the A<r distribution of NLS1 and BLS1 
galaxies (Fig. 3) shows that the two samples are consistent with being drawn from the same 
parent distribution (K-S probability of 0.4). We note that the FWHMs of [SII] and [OIII] core 
correlate well with each other (after removal of 'blue outliers'; Sect. 3.3), demonstrating 
that both emission lines trace similar NLR velocity fields (Fig. 1). Only 40% of the NLS1 
galaxies of our sample have [01] detected. The scatter in crjoi] is relatively small and most 
galaxies are close to the FF05 relation, except two outliers with large api]- 

3.3. Influence of several parameters on the galaxies' location on the M B h — cr 

plane 

We investigated the influence of several quantities - [OIII] CO rc blueshift, L/L^dd, intensity 
ratio FeII4570/H / 3, NLR density, and absolute magnitude - on the galaxies' location on the 
Mbr — c plane. We divided each of these parameters into 3 bins: low values, intermediate 
values, and high values, such that about equal numbers of NLS1 galaxies fell in each bin. We 
measured the velocity shift of [OIII] cor e, ^[oni], relative to [SII] (we use positive velocity values 
for blueshifts). L-^dd was derived from the black hole mass, according to L^dd = 1.3 10 38 
Mbh/M erg/s. The bolometric luminosity Lboi was estimated from the luminosity at 5100A, 
L 510 o, and a bolometric correction of Lboi = 9 AL 510 o- The NLR density was derived from the 
density-sensitive [SII] intensity ratio, [SII]6716/[SII]6731 (Xu et al. 2007), and the absolute 
magnitude from the SDSS i-band magnitude. Our aim is to search for trends across the 
Mbh — o~ plane in order to see whether we can identify physical mechanisms that determine 
an object's location on the plane, and, in particular, that are responsible for the [OIII] 
outliers. 

First, we find that almost all outliers in cqom] actually show a very high blueshift 
of [OIII]core with a velocity shift i>[oni]>150 kms -1 [hereafter referred to as 'blue outliers' 
(Zamanov et al. 2002); marked in Fig. 2]. More generally, we find the trend across the 
M B h — o"[oni] plane that [OIII] lines of higher a also show higher outflow velocity, revealing 
the presence of an extra radial velocity field presumably due to outflows. Among the galaxies 
with positive Act (those below the FF05 relation) 49% are in the highest-velocity bin, while 
only 26% are in the lowest bin. However, if we exclude all galaxies with crromi in the highest 
bin (i>[oiiii>75 kms -1 ), NLS1 and BLS1 galaxies show similar scatter and agree with the 
Mbh — & relation of normal galaxies (Fig. 3; the K-S test shows that the two samples are 
not statistically distinguishable with K-S probability 0.5). We have checked that [SII] is not 
systematically influenced by the outflow components that appear in [OIII] CO rc- NLS1 galaxies 
with blue outliers in [OIII] CO re are uniformly distributed across the Mbh — & plane in [SII]. 
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In order to see whether other parameters affect an object's location on the Mbh — crorni 
plane (Fig. 2), we first removed all objects with the highest [OIII] outflow velocities (f[oini > 
75 km/s), and then studied the correlation of each parameter with the object's deviation 
from the FF05 relation, measured in terms of Act. We cannot identify a strong dependence 
of Act on L/L^&d (Spearman rank correlation coefficient rs=O.lO)0, FeII4570/H/3 (rs=0.05), 
and density (rs=— 0.01). However, if we keep all objects with high blueshifts, LjL^A and 
Act are correlated (rg=0.34). Galaxies with higher absolute magnitudes Mj are, on average, 
located at higher BH masses, as expected, but there is no systematic trend horizontally 
across the plane. 



4. Discussion and conclusions 

We confirm the well-known result (e.g., Grupe 2004) that NLS1 galaxies have, on aver- 
age, lower BH masses than BLS1 galaxies (in our sample, masses range between log -Mbh,nlsi — 
5.7-7.3 M , while logM BH ,BLSi = 6.5-8.4 M ). We further confirm that NLSls are off the 
-Mbh — o~ relation when [OIII] to ti is used as surrogate for a*, and that - after removing blue 
wings - [OIII] core gives more consistent results. However, after correction for blue wings sev- 
eral blue outliers still stick out with very large a values. If these were to reach the M B h — c* 
relation moving vertically on the Mbh — 0" plane by growing their black holes, they would 
have to increase their current BH masses by factors >100, and become more massive than 
the most massive BHs in the BLS1 comparison sample. Instead, the fact that all show high 
outflow velocities in [OIII] cor e demonstrates that the properties of the NLR (unusual outflows 
in the high-ionization NLR) drive the deviation of these NLS1 galaxies from the Mbh — o~ 
relation. Those galaxies which are extreme outliers in [OIII] still follow the Mbh — o" relation 
in [SII]. 

There is good evidence that many NLS1 galaxies are accreting at a high rate relative to 
the Eddington rate (even if not as high as in high-redshift quasars; Warner et al. 2004), so, 
as emphasized by Mathur & Grupe (2005a,b) their BHs must be rapidly growing. Since we 
find that they are on the Mbh — & relation now, their tracks on the Mbh — cr plane must be 
diagonal, else they would deviate in the future. We distinguish two possibilities: (1) NLS1 
galaxies as a class evolve into BLS1 galaxies with respect to their BH mass distribution. 
Then, on average, the BHs of NLS1 galaxies have to grow by a factor ~10. If accreting at 



1 Note that this does not imply an inconsistency with Mathur & Grupe (2005a, b), since their (soft X-ray 
selected) sample includes a significant number of NLS1 galaxies with very high L/L-Edd, above the highest 
value in our sample. 
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L/L^d— 1.0 this would take them 10 8 years. If they did so, a corresponding change in cr* of 
the bulge of the host galaxy is required in order to ensure that the Mbh — o"* relation still 
holds. Recent host galaxy studies of NLSls indicate that they are not preferentially merging 
systems (Krongold et al. 2001, Ryan et al. 2007), but have a relatively high bar fraction 
(Crenshaw et al. 2003, Ohta et al. 2007). In that case, some secular mechanism must be at 
work to adjust the bulge properties. At this stage, we can only speculate that some of the 
feedback mechanisms (Silk & Rees 1998) studied in the context of merger- induced BH fuelling 
(e.g., di Matteo et al. 2005, King 2005) also hold for bar-induced fuelling. Alternatively, 
(2) the NLS1 mass distribution does not evolve into the BLS1 mass distribution. In this 
picture, NLSls are simply low-mass extensions of the BLS1 phenomenon (e.g., McHardy et 
al. 2006). In that case, the high L/L^dd values could just represent a relatively short-lived 
accretion phase which would only insignificantly grow their BHs. Then, within the large 
scatter of a, no adjustment of the host's bulge would be needed. Sources with relatively low 
L/L-Edd could be before or after an accretion event. 

In summary, we conclude that the NLS1 and BLS1 galaxies of our sample, as a class, 
follow the Mbh — o"* relation if we use as surrogate for a* the widths of emission lines which 
are not strongly affected by outflow components ([SII]; or [OIII] core after removing galaxies 
with high outflow velocities v [oni])- Further studies of the evolution of highly accreting AGN 
along the Mbh — o"* relation will shed new light on the evolution of galaxies in general and 
NLSls in particular. 

We thank D. Merritt, H. Zhou, I. Strateva, S. Mathur, D. Grupe, A. Robinson, and 
G. Hasinger for discussions. DX acknowledges the support of the Chinese National Science 
Foundation (NSFC) under grant NSFC-10503005. This research made use of the SDSS data 
base. 
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Fig. 1. — Left panels: Location of NLS1 galaxies (filled symbols) and BLS1 galaxies (open 
symbols) on the M B h — o"[oni] an d -^bh — ^[sii] plane, respectively. In [OIII], black circles 
refer to measurements of the core of the line after correction for asymmetric blue wings, 
while green squares correspond to measurements of the full profile (as described in the text). 
The dashed and solid lines correspond to the M B h — o"* relation of non-active galaxies of 
Tremaine et al. (2002) and of Ferrarese & Ford (2005), respectively, a and FWHM are 
measured in km/s. Right: Correlation between the FWHM of [SII] and [OIII] core (dotted: 
one-to-one line fo equal FWHMs - note the different scales in x and y axis). Blue outliers 
in [OIII] core (f[oni] > 150 km/s) are marked with an extra open square. 
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Fig. 2. — Influence of different parameters on the galaxies' positions on the Mbh — C[oiii] C ore 
plane (open circles: BLS1 galaxies, filled circles: NLS1 galaxies). Each parameter is devided 
into 3 bins and coded by circle size: small: low value, medium: intermediate value, large: 
high value. From left to right: coding according to [OIII] core blueshift (f[oni] <30 km/s, 30- 
75 km/s, >75 kms/s in lowest, intermediate, and highest velocity bin, respectively), L/L Edd 
(<0.6, 0.6-0.9, >0.9), and NLR density (<140 cm" 3 , 140-240 cm" 3 , >240 cm" 3 ). Blue 
outliers in [OIII] corc are marked with an extra open square. 
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Fig. 3. — Histograms of the deviation Aa := log a a b s — log <J pr ed of galaxies from the Mrh — o" 
relation of FF05 (solid line: NLS1 galaxies, dotted line: BLS1 galaxies) for [OIII] CO rc {after 
excluding objects in the highest velocity bin), [OIII] t oti, and [SII], from right to left. 



